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Abstract 

Three-dimensional  computational  simulation  was  employed  to  illustrate  the  performance  characteristics  according  to  the  flow-field  design  by 
solving  the  physics  in  the  flow  field  and  the  diffusion  layer  and  by  calculating  the  electrochemical  reaction  at  the  catalyst  layer.  The  pressure  loss 
and  the  concentration  distribution  in  the  anode  were  analyzed  for  four  types  of  flow  field,  parallel,  serpentine,  parallel  serpentine  and  zigzag  type. 
Also  the  anode  current  density  distribution  was  predicted  at  the  various  overpotentials.  The  cell  performance  was  proportional  to  the  pressure  drop 
for  all  the  flow-field  types.  Zigzag  type  showed  the  best  performance  which  has  a  good  resistance  against  the  fuel  concentration  polarization  and 
the  next  was  serpentine. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Because  of  the  convenience  of  fuel  storage  and  the  com¬ 
pactness,  direct  methanol  fuel  cell  (DMFC)  is  the  most  remark¬ 
able  type  for  potable  applications  among  all  kinds  of  fuel  cells 
[1,2].  DMFC  is  similar  to  proton  exchange  membrane  fuel  cell 
(PEMFC)  in  using  polymer  membrane  as  electrolyte  but  dif¬ 
ferent  in  using  liquid  methanol  as  fuel.  The  fact  that  DMFC 
uses  liquid  fuel  offers  an  advantage  in  fuel  storage  but  brings 
about  other  problems.  One  is  methanol  crossover  to  the  cathode 
side.  Methanol  crossover  is  mainly  occurred  by  diffusion  and 
by  electro-osmotic  drag.  Methanol  crossover  makes  a  mixed 
potential  at  the  cathode  and  therefore  it  decreases  the  cell  per¬ 
formance.  Several  earlier  researchers  have  investigated  about 
methanol  permeation  to  find  out  the  exact  mechanism  and  to 
reduce  methanol  crossover  [3-8].  The  other  problem  is  two- 
phase  flow.  Carbon  dioxide  produced  by  methanol  oxidation  is 
gaseous  phase  although  methanol-water  mixture  is  liquid  phase. 
Therefore,  it  causes  a  big  pressure  loss  in  the  flow  fields  and  also 
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makes  gas  blocks  which  keep  liquid  fuel  from  diffusing  into  the 
catalyst  surface  [9- 1 2] .  In  addition,  it  is  difficult  to  make  uniform 
concentration  because  the  diffusivity  of  methanol  is  smaller  than 
hydrogen  gas.  DMFC  is  greatly  related  not  only  to  a  low  cat¬ 
alytic  activity  of  methanol  electrooxidation  but  to  mass  transfer 
problems.  Anodic,  cathodic  and  overall  reactions  of  DMFC  as 
follows: 


Anode  :  CH3OH  +  H20  -* 

C02  +  6H+  +  6e“ 

(1) 

Cathode  :  |02  +  6H+  +  6e“ 

-  -*  3H20 

(2) 

Overall  :  CH3OH  +  \02  -> 

C02  +  2H20 

(3) 

The  polarization  of  fuel  cell  generally  consists  of  activation, 
ohmic  and  concentration  polarization.  The  shape  of  flow  field 
affects  the  increase  of  the  concentration  overpotential  due  to 
concentration  polarization  which  is  derived  from  the  lack  of 
fuel.  The  effect  of  the  concentration  polarization  is  stronger  at 
the  anode  side  than  at  the  cathode  side  because  of  the  slow  mass 
transfer  of  the  liquid  fuel.  At  the  cathode  side  it  is  more  important 
to  remove  water. 

Kulikovsky  et  al.  [13]  accomplished  two-dimensional  mod¬ 
eling  based  on  mass  and  current  conservation  and  found  out 
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Nomenclature 

c 

molar  concentration  (molm-3  or  M) 

cMeOH 

molar  concentration  of  methanol  (mol  m-3  or  M) 

„MeOH 

c0 

initial  molar  concentration  of  methanol  (mol  m-3 
or  M) 

MeOH 
c  threshold 

molar  concentration  of  methanol  at  threshold 
point  (mol  m3  or  M) 

D 

diffusion  coefficient  (m2  s-1) 

DMeOH.w  diffusion  coefficient  of  methanol  in  water 
(m2  s-1) 

^MeOH  w  effective  diffusion  coefficient  of  methanol  in 
water  (m2  s-1) 

F 

Faraday  constant  (96,487) 

F 

volume  force 

/ 

current  density  (Am-2) 

rMeOH 

70 

exchange  current  density  of  methanol  oxidation 
(Am-2) 

rMeOH 

fe,ref 

exchange  current  density  of  methanol  oxidation 
at  reference  point  (A  m-2) 

k 

permeability  (m2) 

n 

a  normal  vector  to  the  boundary 

^reaction 

flux  derived  by  reaction  (molm-2  s-1) 

P 

pressure  (Pa) 

P  interface 

pressure  at  the  interface  between  channel  and 
backing  layer  (Pa) 

R 

gas  constant  (8.3 14  J  mol-1  K-1) 

t 

time  (s) 

T 

temperature  (K) 

u 

velocity  vector  (ms-1) 

Uo 

initial  line  velocity  (ms-1) 

Greek  letters 

O'a 

anodic  transfer  coefficient 

£ 

porosity 

0* 

anodic  overpotential  (V) 

Ji 

viscosity  (kgm-1  s-1) 

P 

density  (kgm-3) 

that  the  transport  of  methanol  was  determined  by  the  pressure 
gradient  and  that  the  creation  of  ‘Shaded’  zone  was  dependent 
on  the  hydraulic  permeability  of  the  backing  layer.  Wang  and 
Wang  [14]  simulated  two-phase  flow  in  liquid-feed  DMFC  by 
the  two-dimensional  model  which  included  the  calculation  of 
the  current  density  distribution  and  the  methanol  crossover  flux 
and  the  prediction  of  the  unit  cell  performance.  Scott  and  Argy- 
ropoulos  [15]  studied  one-dimensional  model  of  DMFC  focused 
on  electrocatalysis  to  predict  anode  polarization  behavior.  Dohle 
and  Wippermann  [16]  investigated  the  anode  overpotential,  the 
cathode  overpotential  and  the  methanol  permeation  in  order 
to  experimentally  obtain  an  empirical  equation.  Their  model 
described  the  current  density  characteristics  and  the  methanol 
permeation.  Although  the  above  models  were  mainly  mathe¬ 
matical  approaches  about  the  electrochemical  reaction  and  the 
mass  transfer,  they  did  not  include  the  effect  of  flow  fields.  On  the 


researches  about  the  flow-field  pattern,  Arico  et  al.  [17]  investi¬ 
gated  the  features  of  the  serpentine  and  interdigitated  flow-field 
type  for  DMFC.  According  to  their  experiments,  the  interdigi¬ 
tated  type  significantly  enhanced  the  mass  transport  and  hence 
showed  the  higher  maximum  power  outputs  compared  to  the 
classical  serpentine  geometry.  Tiiber  et  al.  [18]  measured  the 
performance  of  PEMFC  and  DMFC  for  serpentine,  parallel  and 
fractal  type  of  flow  field  as  the  flow-field  design  and  the  ser¬ 
pentine  type  showed  the  best  result.  They  speculated  that  the 
serpentine  channel  forces  the  products  out  of  the  cell  by  its  higher 
pressure  drop  and  that  the  parallel  and  fractal  designs  can  suffer 
from  the  blocked  channels  by  inhomogeneous  flow  distribution. 

In  this  work  three-dimensional  computational  simulation  was 
used  in  order  to  illustrate  the  causes  of  performance  difference 
according  to  the  flow-filed  design.  Parallel,  serpentine,  paral¬ 
lel  serpentine  and  zigzag  shapes  were  simulated  and  estimated 
for  the  flow  field  of  the  anode  by  comparing  the  concentration 
distribution  and  the  current  density  distribution  at  the  various 
overpotentials. 

2.  Experimental 

FEMLAB,  a  commercial  program,  is  used  in  this  work. 

2.7.  Assumptions 

1 .  The  DMFC  system  is  steady  state. 

2.  Gravity  is  ignored. 

3.  The  system  is  isothermal. 

4.  Two-phase  flow  due  to  carbon  dioxide  is  ignored. 

5.  Vaporization  and  crossover  of  the  anode  fuel  are  ignored. 

6.  The  flow  is  laminar  and  incompressible. 

7.  Fuel  consumption  does  not  affect  Navier-Stokes  momentum 
balance  in  the  fuel  channel. 

8.  Thickness  of  the  catalyst  layer  is  zero. 

9.  Fuel  concentration  change  due  to  water  consumption  is 
ignored. 

2.2.  Geometry 

Fig.  1  shows  the  components  of  DMFC.  Flow-field  plate  plays 
the  part  of  the  fuel  flow  field  and  of  the  current  collector  from 
membrane  electrode  assembly  (ME A).  Fuel  flows  through  the 
fuel  channels  and  diffuses  into  the  backing  layer.  The  backing 
layer  supports  the  catalyst  layer  and  controls  the  fuel  diffusion 
into  the  catalyst  layer.  Electrochemical  oxidation  and  reduction 
occurred  at  the  catalyst  layer  and  membrane  plays  a  part  in  the 
proton  conduction.  The  simulated  part  in  this  study  is  the  anode 
side,  the  inside  of  the  red  box  in  Fig.  1.  An  overview  of  the 
geometry  is  shown  in  Fig.  2. 

Four  types  of  the  flow  field  which  is  parallel,  serpentine,  par¬ 
allel  serpentine  and  zigzag  were  used  for  the  simulation.  In  all 
cases,  fuel  input  and  fuel  output  are  located  diagonal  to  each 
other  and  the  electrode  area  is  about  2.25  cm2.  Height  and  width 
of  the  channel  are  fixed  as  1  mm.  Fuel  flows  according  to  the  flow 
fields  are  described  in  Fig.  3  and  the  geometrical  specifications 
are  shown  in  Table  1 . 
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Simulation  Part  - 

Fuel  Channel 

Anode  Catalyst  Layer 
Membrane 
Cathode  Catalyst  Layer 

Air/02  Channel 


Fig.  1.  Cross-sectional  view  of  a  DMFC  and  the  simulation  part. 


Fig.  2.  An  overview  of  the  geometry  (serpentine  type). 


Table  1 


Specifications  of  the  flow  fields 


Electrode 
area  (cm2) 

Width  x 

length 

(cm2) 

Channel 

area 

(cm2) 

Rib 

area 

(cm2) 

Rib 

area 

ratio 

Parallel 

2.25 

1.5  x  1.5 

1.41 

0.84 

0.373 

Serpentine 

2.21 

1.7  x  1.3 

1.25 

0.96 

0.434 

Parallel  serpentine 

2.255 

2.05  x  1.1 

1.32 

0.935 

0.415 

Zigzag 

2.2575 

1.5  x  1.55 

1.3875 

0.87 

0.385 

2.3.  Governing  equations 

In  the  flow-field  channel,  Navier-Stokes  equations  (Eq.  (4)) 
were  employed  to  calculate  velocity  and  pressure.  The  momen¬ 
tum  balance  and  the  continuity  equation  are  as  follows: 

p - V  •  /t(Vu  +  (Vu)T)  +  p(u  •  V)u  +  VP  =  F  (4) 

dt 

V  •  u  =  0  (5) 

where  F  denotes  the  volume  force  term  such  as  gravity  force. 

On  the  other  hand,  the  momentum  balance  in  the  porous  back¬ 
ing  layer  can  be  described  by  Darcy’s  law, 

k 

u  =  -- VP  (6) 

The  convective  velocity  is  determined  by  pressure  gradient 
because  the  permeability  and  the  viscosity  are  constant.  Pressure 
gradient  can  be  derived  from  Navier-Stokes  equations  because 
the  flow-field  channel  and  the  backing  layer  are  contacted  each 
other. 


Fig.  3.  Fuel  flows  according  to  the  flow-held  type:  (a)  parallel,  (b)  serpentine,  (c)  parallel  serpentine  and  (d)  zigzag. 


878 


M.-s.  Hyun  et  al.  /  Journal  of  Power  Sources  157  (2006)  875-885 


Table  2 


Physical  and  chemical  properties 


Parameter 

Value 

Reference 

Anodic  transfer  coefficient, 

0.239 

[14] 

Oia 

Reference  exchange  current 

94.25 

[14] 

density  of  methanol 
oxidation  at  80  °C, 

C?H  (Am_2) 

Cell  temperature,  T  (K) 

353 

Channel  height,  h  (m) 

0.001 

Channel  width,  w  (m) 

0.001 

Density  of  1  M  methanol,  p 

980 

Measured 

(kgm-3) 

Diffusion  coefficient  of 

jQ— 5.4163  —  999.778/7" 

[14] 

methanol  in  water, 

DMeOH  (m2  s-1) 

Input  methanol 

1000 

concentration,  c^e0H 
(molm-3) 

Input  velocity  of  fuel,  v 

1.5  X  10“8 

(m3  s-1) 

Permeability  of  backing 

i  x  nr11 

Assumed 

layer,  k  (m2) 

Porosity  of  backing  layer,  s 

0.77 

Assumed 

Thickness  of  backing  layer, 

0.0003 

Assumed 

l  (m) 

Viscosity  of 

0.458509-5.304741 

[19] 

methanol-water  mixture, 

x  10“37’+  2.31231 

U  (kgm-1  s-1) 

x  10“5r2  -4.49161  x  io~8r3 
+  3.27681  x  KT11:?4 

The  convection  and  diffusion  equation  is, 
dc 

f  =  V-(DVc-cu)  (7) 

ot 

The  expression  within  the  parenthesis  represents  the  flux,  where 
the  first  term  describes  the  transport  by  diffusion  and  the  second 
represents  the  convective  flux.  The  velocity  of  the  convective  flux 
was  obtained  by  Navier-Stokes  equations  or  Darcy’s  equation. 

The  bottom  of  the  backing  layer  is  regarded  as  the  catalyst 
layer.  In  the  catalysts  layer  only  the  electrochemical  kinetics  are 
considered  because  the  thickness  of  the  catalyst  layer  is  assumed 
as  zero.  In  this  work,  Tafel  kinetic  equation  of  methanol  oxida¬ 
tion  was  used  which  was  proposed  by  Wang  and  Wang  [14]. 
According  to  Tafel  equation,  methanol  oxidation  is  a  zero-order 
reaction  when  the  methanol  concentration  is  higher  than  0. 1  M. 
But  under  the  threshold  value,  0. 1  M,  the  methanol  oxidation 
is  considered  as  a  first-order  reaction  and  the  exchange  current 
density  decreases  rapidly  (Table  2). 


rMeOH 

70 


exp 


rMeOH  _  rMeOH 
70  —  70,ref 


cMeOH  \n 

MeOH  ] 

L  threshold  / 


n  = 


0; 

1; 


/„MeOH  MeOH  ^ 
'T  —  c  threshold' 

/  JMeOH  .  MeOH  x 
'T  c  threshold' 


MeOH 

^threshold 


0.1  M 


(10) 


2.4.  Boundary  conditions  and  parameters 

2.4.1.  Navier-Stokes  equations 
At  wall, 

u  =  0  (11) 

At  the  entrance  of  fuel  input, 

u  •  n  =  uo  (12) 

where  uo  is  the  line  velocity  of  fuel. 

At  the  exit  of  fuel  output, 

P  =  0  (13) 


2.4.2.  Darcy's  law 
At  wall, 


k 

/a 


V  p  •  n  =  0 


(14) 


At  the  interface  between  the  flow-field  channel  and  the  backing 
layer, 


P  —  ^interface  (13) 

where  interface  Is  derived  from  Navier-Stokes  equations. 

2.4.3.  Convection  and  diffusion  equation 

The  diffusion  coefficient  in  the  porous  backing  layer  is  cor¬ 
rected  with  the  porosity  of  the  backing  layer  by  Bruggemann’s 
correction, 

^MeOH,w  =  ^MeOH,w  X  (16) 

The  effective  diffusion  coefficient  was  used  in  the  backing  layer 
domain  and  the  bulk  diffusion  coefficient  was  used  in  the  fuel 
channel  domain. 

The  velocity  from  Navier-Stokes  equations  was  used  in  the 
channel  domain  and  the  convection  velocity  from  Darcy’s  law 
was  used  in  the  backing  layer. 


At  wall, 


(D  Vc  —  cu)  •  n  =  0 


(17) 


At  the  entrance  of  the  fuel  input, 


c  =  c 


MeOH 

0 


where  c^eOH  is  the  initial  methanol  concentration. 
At  the  exit  of  the  fuel  output, 


(18) 


n  •  (DMeOH,WVc)  =  0  (19) 

At  the  bottom  of  the  backing  layer, 

(^MeOH,wVc  “  CU)  ’  H  —  Nreaction  (20) 

where  Nreaction  is  the  methanol  flux  which  is  obtained  from  the 
methanol  oxidation  kinetics. 
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2.4.4.  Kinetic  equation  of  methanol  oxidation 

The  exchange  current  density  (9)  can  be  rewritten  as  follows: 


/  iivieun  /  ivieun  \/ 

rMeOH  _  rMeOH  ,  V0,ref  /  c  threshold  A1 L 

Q  *  0  rpf  * 


MeOH /JVIeOH  \/ JVIeOH  JVIeOH  ^  i  /  zMeOH  /  JVIeOH  v  JVIeOH MeOH 


—  c 


threshold 


\  _  |  (  nvieun  / 

1  I 'To, ref  /c 


threshold 


)(cJ 


threshold 


)l 


0,ref 


Therefore,  Tafel  kinetic  equation  (Eq.  (8))  can  be  rewritten  as 
follows: 


MeOH  /  JVIeOH  \/ JVIeOH  JVIeOH 


/rivieun  /  ivieun 

T  _  |  rMeOH  i  'To, ref  /  cthreshold^c 

i  —  I  in  i 


—  C 


MeOH  /JVIeOH  \/ JVIeOH  ^.MeOH 


threshold 


x  _  i/  rivieun  / 

1  I  Do, ref  !C 


threshold 


)(cJ 


—  C 


threshold 


)l 


«a  F 


0,ref 


exp 


RT 


?7a 


(21) 


(22) 


and  the  methanol  flux  is 

1 


rMeOH  /  JVIeOH  w  JHeOH 


N 


reaction  — 


6F 


rMeOH  I 
70,ref 


(  nvieun  / 
t^O^ef  /cthreshold 


)(cJ 


MeOH  \  _  |  (  fivieun  /  ivieun  \/  ivteuri  _  ivieun  \i 
c threshold^  It^Q^ef  ' c threshold cthreshold^l 


rMeOH  /  ^MeOH 


MeOH  ^MeOH 


exp 


a^F 


RT 


(23) 


3.  Result  and  discussions 

The  calculated  methanol  concentration  distribution  in  the 
overall  geometry  of  the  serpentine  type  is  shown  in  Fig.  4  when 
the  overpotential  of  0.45  V  is  imposed  to  the  electrode.  The 
concentration  difference  between  the  flow  field  and  the  cata¬ 
lyst  layer  is  big  while  the  difference  between  fuel  input  and 
output  in  the  flow  fields  is  not  big.  Particularly,  concentration  at 
rib  areas  of  the  catalyst  layer  which  contact  with  bipolar  plate 
is  very  low  because  fuel  cannot  contact  with  the  backing  layer 
directly.  Methanol  concentration  at  the  catalyst  layer  is  much 
lower  than  in  the  flow  fields  and  there  is  also  a  wide  difference 
of  methanol  concentration  among  the  positions  on  the  catalyst 
layer. 

When  the  overpotential  of  0.3  V  is  imposed,  the  methanol 
concentration  distributions  at  the  catalyst  layers  are  shown  in 
Fig.  5.  The  minimum  methanol  concentrations  are  higher  than 
100  mol  m-3,  the  threshold  value,  for  all  the  flow-field  types, 
although  low  concentrations  which  are  expressed  in  blue  in  Fig.  5 


are  shown  at  ribs  for  parallel  and  parallel  serpentine  type.  There¬ 
fore,  there  are  no  current  deviations  within  the  all  catalyst  layers 
and  the  calculated  current  densities  for  all  the  flow-field  types 
are  995  A  m-2  equally. 

Fig.  6  shows  the  methanol  concentration  distributions  when 
the  overpotential  of  0.4  V  is  imposed.  The  zigzag  type  shows  the 
highest  average  concentration  and  the  parallel  type  is  the  low¬ 
est.  Comparing  Fig.  5,  overall  methanol  concentration  decreased 
greatly  due  to  the  increment  of  oxidation  kinetics.  Methanol 
concentration  of  parallel  type  reaches  to  the  threshold  value 
on  the  whole  areas  of  the  catalyst  layer.  Fig.  7  shows  the 
current  density  distributions  when  the  overpotential  of  0.4  V 
is  imposed.  Dead  zone  is  observed  at  the  ribs  of  the  par¬ 
allel  and  parallel  serpentine  types  partially  although  the  ser¬ 
pentine  and  zigzag  did  not  yet  fall  into  the  concentration 
polarization. 

In  the  case  of  0.45  V  overpotential  (Figs.  8  and  9),  most 
sections  of  the  parallel  and  the  parallel  serpentine  reach  to  the 
threshold.  The  serpentine  type  starts  to  fall  into  the  concentration 


Slice:  Concentration,  c 


Max:  1000 
1000 

900 

800 
700 

600 
500 
400 
300 
200 
100 
0 

Min:  0 


Fig.  4.  The  distribution  of  the  calculated  methanol  concentration  for  the  serpentine  type  when  the  overpotential  of  0.45  V  is  imposed  (unit:  molm  3). 
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Fig.  5.  Concentration  distributions  at  the  electrode  surface  when  the  overpotential  of  0.3  V  is  imposed.  Average  concentrations  on  the  catalyst  layer:  (a)  525  molm 
for  parallel,  (b)  780 molm-3  for  serpentine,  (c)  676 molm-3  for  parallel  serpentine  and  (d)  803  molm-3  for  zigzag. 


polarization  although  zigzag  type  is  still  keeping  the  theoretical 
current  density  without  concentration  polarization.  When  0.5  V 
overpotential  is  imposed  (Figs.  10  and  11),  a  significant  lack 
of  methanol  appeared  in  all  types  of  the  flow  field  and  also  the 


difference  in  the  current  densities  increased  within  each  catalyst 
layer.  The  concentration  polarization  generates  easily  at  ribs  and 
this  area  acts  as  ‘dead  zone’  at  which  current  density  decreases 
rapidly.  In  case  of  the  zigzag  and  the  serpentine  types,  relatively 


Fig.  6.  Concentration  distributions  at  the  electrode  surface  when  the  overpotential  of  0.4  V  is  imposed.  Average  concentrations  on  the  catalyst  layer:  (a)  178  mol  m 
for  parallel,  (b)  517 molm-3  for  serpentine,  (c)  338  molm-3  for  parallel  serpentine  and  (d)  570 molm-3  for  zigzag. 
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Fig.  7.  Current  distributions  at  the  electrode  surface  when  the  overpotential  of  0.4  V  is  imposed.  Average  current  densities  on  the  catalyst  layer:  (a)  1827  Am  2  for 
parallel,  (b)  2184  Am-2  for  serpentine,  (c)  2086  Am-2  for  parallel  serpentine  and  (d)  2184  Am-2  for  zigzag. 


uniform  concentration  distribution  was  found  on  the  catalyst 
layers. 

The  velocity  profile  at  the  bottom  of  the  backing  layer  shows 
the  reason  why  the  serpentine  type  shows  a  relatively  uniform 


fuel  concentration.  Fig.  12  shows  the  velocity  distributions  at 
the  bottom  of  the  backing  layer  of  the  serpentine  and  the  parallel 
type  flow  channel.  The  velocity  of  parallel  type  is  of  the  order 
of  10~7ms_1  which  means  the  convective  velocity  is  nearly 


Fig.  8.  Concentration  distributions  at  the  electrode  surface  when  the  overpotential  of  0.45  V  is  imposed.  Average  concentrations  on  the  catalyst  layer:  (a)  84  mol  m 
for  parallel,  (b)  303  molm-3  for  serpentine,  (c)  166  mol  m~3  for  parallel  serpentine  and  (d)  363  molm-3  for  zigzag. 
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Fig.  9.  Current  distributions  at  the  electrode  surface  when  the  overpotential  of  0.45  V  is  imposed.  Average  current  densities  on  the  catalyst  layer:  (a)  2087  Am-2  for 
parallel,  (b)  3152  Am-2  for  serpentine,  (c)  2723  Am-2  for  parallel  serpentine  and  (d)  3235  Am-2  for  zigzag. 


negligible  and  therefore  the  diffusion  through  the  backing  layer 
is  mainly  controlled  by  the  mass  transfer  in  the  backing  layer. 
However,  in  serpentine  type  the  velocity  at  the  bottom  of  the 
backing  layer  is  103  or  104  times  higher  than  that  in  the  paral¬ 


lel  type.  It  means  that  both  convection  and  diffusion  control  the 
mass  transfer  in  serpentine  type.  The  convective  velocity  is  orig¬ 
inated  from  the  pressure  gradient  of  the  flow  fields  and  makes 
the  mass  transfer  of  fuel  easy.  Therefore,  the  velocity  from  the 
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Fig.  10.  Concentration  distributions  at  the  electrode  surface  when  the  overpotential  of  0.5  V  is  imposed.  Average  concentrations  on  the  catalyst  layer:  (a)  47  mol  m 
for  parallel,  (b)  107  mol  m-3  for  serpentine,  (c)  66  mol  m-3  for  parallel  serpentine  and  (d)  119  mol  m-3  for  zigzag. 
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Fig.  1 1 .  Current  distributions  at  the  electrode  surface  when  the  overpotential  of  0.5  V  is  imposed.  Average  current  densities  on  the  catalyst  layer:  (a)  2186  A  m  2  for 
parallel,  (b)  3891  Am-2  for  serpentine,  (c)  3003  Am-2  for  parallel  serpentine  and  (d)  4283  Am-2  for  zigzag. 


Fig.  12.  Velocity  profile  at  the  bottom  of  the  backing  layer:  (a)  serpentine  and 
(b)  parallel  (unit:  ms-1). 


position  A  to  the  position  B  in  Fig.  12(a)  is  bigger  than  that 
from  the  position  C  to  the  position  D  in  Fig.  12(a).  Zigzag  type 
compensates  this  kind  of  the  velocity  difference.  Finally  high 
pressure  drop  contributes  the  formation  of  uniform  concentra¬ 
tion.  The  parallel  type  does  not  need  high  pressure  drop  in  the 
channel  however  the  serpentine  and  the  zigzag  type  require  the 
high  pressure  drop.  The  pressure  difference  between  the  parallel 
and  the  serpentine  types  reaches  about  102  times  (Fig.  13). 

The  parallel  serpentine  type  is  an  intermediate  type  between 
the  serpentine  and  the  parallel.  The  channel  shape  of  the  parallel 
serpentine  is  meandering  which  is  similar  to  serpentine  type  and 
the  parallel  channel  of  the  parallel  serpentine  is  similar  to  the 
parallel  type.  The  zigzag  type  is  based  on  the  serpentine  type. 
The  main  idea  of  the  zigzag  type  is  to  make  the  width  from  C  to 
D  close  because  the  fuel  concentration  near  the  bended  point  of 
channel  is  low. 

Fig.  14  presents  an  anode  overpotential  curve  according  to  the 
flow-field  design.  Each  point  is  originated  from  the  average  cur¬ 
rent  densities  calculated  from  the  methanol  concentration  at  the 
catalyst  layer  when  the  corresponding  overpotential  is  imposed. 
Dot  line  is  a  theoretical  activation  polarization  curve  in  which  the 
concentration  polarization  is  not  considered.  According  to  the 
simulated  results,  all  types  of  the  flow  fields  agree  the  theoretical 
values  well  until  the  overpotential  of  0.35  V  whereas  they  can¬ 
not  overcome  the  concentration  polarization  and  finally  bring 
out  the  sudden  increment  of  overpotential.  The  overpotential 
curves  of  the  parallel,  the  parallel  serpentine,  the  serpentine  and 
the  zigzag  types  deviate  the  theoretical  overpotential  curve  from 
142  mA  cm-2  at  0.35  V,  208  mA  cm-2  at  0.4  V,  3 15  mA  cm-2  at 
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Fig.  13.  Pressure  drop  in  the  flow-field  channel  when  the  fuel  input  velocity  is  1.5  x  10  8  m3  s  1 :  (a)  0. 18  Pa  for  parallel,  (b)  19.70  Pa  for  serpentine,  (c)  4.95  Pa  for 
parallel  serpentine  and  (d)  22.50  Pa  for  zigzag. 


Fig.  14.  Predicted  anode  overpotential  curve  according  to  the  flow-field  design. 

0.45  V  and  391  mAcm~2  at  0.475  V,  respectively.  Zigzag  type 
can  be  estimated  as  the  best  flow-field  shape  in  this  work. 

4.  Conclusions 

A  three-dimensional  computational  simulation  of  DMFC 
anode  was  developed  for  the  parallel,  the  serpentine,  the  par¬ 
allel  serpentine  and  the  zigzag  types  in  order  to  illustrate  the 
performance  difference  according  to  the  flow-field  design.  By 
examining  the  methanol  concentration  distribution  and  the  cur¬ 
rent  density  distribution  the  performance  of  the  anode  for  each 


flow  field  was  estimated.  The  zigzag  type  showed  the  best  perfor¬ 
mance  which  has  a  good  resistance  against  the  fuel  concentration 
polarization  and  the  next  are  the  serpentine,  the  parallel  ser¬ 
pentine  and  the  parallel  in  order.  The  cell  performance  was 
proportional  to  the  pressure  drop  for  all  the  flow-field  types. 
Pressure  energy  which  forces  liquid  fuel  out  of  the  channels 
helps  fuel  to  spread  uniformly  by  driving  the  convective  veloc¬ 
ity  in  the  backing  layer.  In  this  work  zigzag  type  showed  the  best 
anode  performance  whereas  needed  the  highest  pressure  drop. 
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